Nanoscale ferromagnet-superconductor-ferromagnet switches controlled by 

magnetization orientation 
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We study clean ferromagnet-superconductor-ferromagnet (FSF) nanostructures in which the mag- 
netization of the F layers can be parallel (P) or antiparallel (AP). We consider the case where the 
thickness of the S layer is of order of the coherence length, with thinner F layers. We find that 
reversing the direction of the magnetization in one of the F layers leads in general to drastic changes 
in the superconductor's state. Under a wide variety of conditions, the AP geometry favors supercon- 
ductivity. Magnetization reversal in one of the F layers can lead to the superconductivity turning 
on and off, or to switching between different states. Our results are obtained via self consistent 
solution of the Bogoliubov-de Gennes equations and evaluation of the condensation energies of the 
system. 

PACS numbers: 74.45.+C, 74.25.Fy, 74.78.Fk 



Introduction: Within the emerging field of 
spintronicsi considerable interest has developed in 
devices in which proximity effects are used to control 
the superconductivity via the spin degree of freedom 
in ferromagnet (F) and superconductor (S) layered 
systems. A large part of the motivation for this interest 
follows from earlier studies of systems that involve non- 
magnetic normal metal layers sandwiched between two 
ferromagnetic layers^ (FNF geometry). In such devices 
the resistance of the system can change substantially in 
the presence of a perturbing magnetic field. This change 
mainly arises from the spin-dependent scattering at the 
interfaces. The ensuing giant magnetoresistance (GMR) 
effect is found in spin-valves and magnetic multilayers 
where the relative orientations of the magnetization in 
alternate ferromagnetic layers change as a function of 
an applied field. If the local magnetization orientations 
are antiparallel (AP) the scattering will be stronger for 
a particular spin component, but if the magnetization 
vectors are aligned the more weakly scattered spin 
component carries the current with a lower resistivity. 

If the nonmagnetic insert is replaced by a thin super- 
conductor, resulting in a ferromagnet-superconductor- 
ferromagnet (FSF) junction, a different type of spin- 
valve or spin-switch can be created^. The proximity ef- 
fects arising from the mutual influence of the magnetic 
and superconducting order parameters embody a variety 
of novel spin-valve effects and device concepts, includ- 
ing high density nonvolatile memory, 4 and magnetic sen- 
sors. The mechanism behind such devices is ultimately 
based^ on the damped oscillatory nature of the super- 
conductor order parameter in the F regions, and the 
associated magnetic correlations and destruction of su- 
perconductivity in the S layer. In the transport regime, 
and with AP alignment of the magnetizations in the F 
layers, a nonequilibrium spin density can accumulate in 
the superconductonSi^, destroying the gap and result- 



ing in a higher resistance state for a given temperature £ 
Thus the superconducting correlations are controlled by 
the relative magnetization orientation in the F layers. 
Also, quasiclassical thermodynamic considerations indi- 
cate that the transition temperature, T c , can be modified 
in a controlled way, thus allowing supercurrent to flow in 
a predictable manner j^iiSiiiii^ yielding another type of 
spin switch. The superconducting order parameter and 
T c in this case are again greatest when the the magnets 
are in the AP configuration, a result shown to hold for 
atomic thickness FSF layers as welli^ When the super- 
conducting system goes normal, the phase transition is 
second order for AP magnetizations in the F layers and 
can be first order for parallel (P) magnetizations if the 
F layers are thin enough and the interface transparency 
is high. 14 If the outer ferromagnets are semiconducting 
insulators, the T c variations have different signatures de- 
pending on whether the superconductor is in the singlet 
or triplet stated An absolute spin- valve effect can occur 
at spin-active interfaces in which the tunneling current is 
finite for a range of voltages^ 

These types of devices are in general most effective, 
and the effects most pronounced, for junctions with 
clean interfaces and thin superconductors, 11 The litho- 
graphic, sputtering, and epitaxial methods used in spin- 
switch fabrication permit the creation of structures as 
thin as a few atomic layers that have atomically flat in- 
terfaces. Moreover, high quality magnetic and nonmag- 
netic metallic films with an electron mean free path ex- 
ceeding 150 A can also be readily fabricated. One of 
the earliest experiments using FSF junctions involved 
CuNi/Nb/CuNi films, and a magnetization direction de- 
pendence on T c was reported £ A spin switch was re- 
cently investigated using Lao.7Cao.3Mn03/YBa2Cu307 
superlattices that had large magnetoresistance when in 
the superconducting stated Spin valve core structures 
involving Nb/CuNi sandwiches have also very recentljii& 
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been reported. It is possible, in an FSF sandwich with 
AP magnetizations, for the electron in one of the magnets 
to be Andreev reflected as a hole of the opposite spin in 
the other ferromagnetic This process of crossed Andreev 
reflection is believed to be behind the results of exper- 
iments involving subgap transport in Al/Fe hybrids. 18 
An enhancement of the critical current and T c in Nb/Co 
was observed and was attributed to a reduced exchange 
interaction in the domain structure of Coi^ Likewise, 
a type of spin-switch involving Nb/Pcrmalloy layers re- 
vealed through transport measurements, a decrease in 
the suppression of superconductivity^ It was argued 
that the superconductivity is increased when the mag- 
netic domains are oriented differently, effectively aver- 
aging in a way that reduces the effects of the exchange 
field. 

Spurred by these important advances, we investigate 
here the effect of reversing one of the magnetizations in 
clean FSF nanojunctions. We consider the relevant case 
where the coupling between the S and F regions is ap- 
preciable, namely a thin S layer (of order of the BCS 
coherence length £o) and relatively large magnetic ex- 
change fields. Our results are based upon numerical self 
consistent solution of the microscopic Bogoliubov de- 
Gennes (BdG) equations. This method is most appropri- 
ate for the situation described above. We calculate the 
pair potential A(r), the condensation energy, and the 
local density of states (LDOS) for both the P and AP 
magnetization configurations, over a range of values of 
the relevant parameters. Our analysis will demonstrate 
that under many conditions the system can be made to 
switch from a superconducting state to a normal one, 
at low temperatures, by flipping the (collinear) magne- 
tization orientation in one of the F layers, which can be 
achieved via an applied field. This will be illustrated by 
calculating the condensation energy as a function of ferro- 
magnet thickness and exchange energy. We find that the 
AP state is always the lowest energy state, and thus the 
most favorable. We conclude that the proximity effects 
that occur with P magnetization in successive F layers 
become substantially modified when adjacent F layers 
have AP magnetization alignment. The pair amplitude 
and LDOS also display experimentally discernible char- 
acteristics that depend on whether the magnets are in 
the P or AP configuration. 

Methods: The equations relevant to the microscopic 
theory of inhomogeneous superconductivity are the 
Bogoliubov-de Gennes (BdG) equations^ We consider 
here an FSF structure translationally invariant in the 
x — y plane, with interfaces normal to the z direction. 
We assume parabolic bands with bandwidths Ep in the 
S layer and Ep± = Ep±ho in the magnet, where ho is the 
Stoner exchange field. In the P geometry the sign of ho is 
the same in both layers, while in the AP geometry it is 2 ^ 
the opposite. The dimensionless parameter / = ho/Ep 
characterizes the magnetic strength. We include interface 
scattering characterized by delta functions of strength H 
(dimensionless strength Hb = mH/kp). We have writ- 
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FIG. 1: (Color online) The spatially dependent pair ampli- 
tude (normalized to its T = bulk value) for a FSF trilayer, 
plotted vs. Z = Ufz. The contrasting cases where the magne- 
tizations in the F layers are parallel (label P), or antiparallel 
(label AP) are shown. Because of symmetry, only half of the 
sample is included. Each magnet has width Df = kpdLF = 10. 
The dashed vertical line is at the F/S interface. Two values 
of the scattering strength Hb are considered. For the smaller 
value the pair amplitude vanishes in the P configuration. The 
pair amplitude is always larger in the AP configuration. The 
exchange parameter is / = 0.2. 

ten the BdG equations in this geometry and with these 
assumptions in previou o 21 i 24 ! work, where we have also 
discussed^ 5 - the specific numerical methodology we use. It 
is therefore unnecessary to repeat these derivations here. 
The exact quasiparticle energies and amplitudes are thus 
obtained by repeated iteration of the BdG equations and 
the associated self consistency condition for A(z). 

Once the energy spectra and pair potential are found, 
the condensation free energy, (or, in the T —> limit 
the condensation energy) can be calculated. This is in 
principle straightforward, although numerically quite dif- 
ficult. We use2£ a particularly convenient approachS&Si 
which yields for the condensation energy AEq the result: 

A^E^-E^+r^^' a) 

where g is the BCS coupling constant in S, d the to- 
tal sample width, e„ and e°; are the free-particle energy 
spectra corresponding respectively to the superconduct- 
ing and normal (A(z) = 0) states, and the indices de- 
note the appropriate 25 quantum numbers. The sums are 
performed over energies less than the usual cutoff ujd- 
Similarly, from the calculated self-consistent eigenvalues 
and eigenfunctions one can calculate the LDOS, N(z, e). 
This quantity is discussed below. 

Results: In the geometry we consider, the inner su- 
perconductor layer of width ds is sandwiched between 
two ferromagnet layers of equal width, dp. The thick- 
ness ds is chosen to be ds = £o: and we take kp^o = 100 
and u>d — 0.04-Ep. All results correspond to low temper- 
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FIG. 2: (Color online) The condensation energy, AEq, nor- 
malized to 7V(0)Aq, plotted vs / for two values of Hb- The 
oscillatory behavior is complicated, but for both Hb values 
the AP state is favored over the whole I range. 



ature, T = 0.01T C . 

In Fig. we plot the pair amplitude (the average 

( - 0x(r)V-'T( r )): where the are the usual annihilation op- 
erators), normalized to its zero temperature bulk value, 
as a function of dimensionless distance. Results are plot- 
ted both for both the P and AP cases. Two values of the 
interface scattering parameter are considered, a small one 
{Hb = 0.15) when the proximity effect is strong, and a 
larger one (Hb = 0.6) when it is weaker. It is clear that 
the results depend crucially on the relative magnetization 
orientation, with the superconducting amplitude being 
weakest in the P case. The effect is magnified for the 
smaller Hb value, where interface scattering is reduced. 
In that case one can see in the figure that the pair am- 
plitude vanishes in the P case, while being quite robust 
in the AP situation. The superconductivity can thus be 
switched on and offhy reversing the magnetization in one 
of the F layers. 

This favoring of the AP configuration is, qualitatively, 
a very general phenomenon. To see this, it is very conve- 
nient to analyze the pair condensation energy (Eq. Q). 
The trends in the pair amplitude, such as those in Fig.^ 
should be reflected in the condensation energy, which 
should then be lower (higher in absolute value) in the AP 
case than in the P configuration. AEo is shown in the 
next two figures, normalized to N(0)Aq which is twice 
its bulk zero temperature value. In Fig. |2 this normal- 
ized quantity is plotted as a function of / for two values 
of Hb- The F width is kept fixed at Dp = kpdp = 10 
(recall that D$ = kp^Q always). The entire range of I 
from the nonmagnetic (I = 0) limit to the half metal- 
lic case (I = 1) is spanned. For all nonzero I, the AP 
case is always favored. This trend persists even for larger 
Hb (not shown) where the proximity effect is weaker. At 
1=1 only one spin band is populated at the Fermi level 
and consequently |A_E | is large, as Andreev reflection is 
depressed and the Cooper pairs are more restricted to the 




FIG. 3: (Color online) The normalized condensation energy 
vs. the dimensionless ferromagnet width Df- Two values of / 
are considered, as indicated, at Hb = 0.3. For a given I, the 
AP state has lower AEq over the whole Df range included. 

S region. We also see in the figure that the difference in 
condensation energies for P and AP geometries at fixed 
thickness is a weak function of I, except at small /. 

The AP configuration continues to be preferred when 
the thickness Dp is varied at constant I. This is shown 
in Fig. |3 where we plot the normalized AEq versus Dp. 
Results for two values of I are plotted, and both P and 
AP configurations are studied. As Dp —> 0, one is left 
only with the superconductor and AEq approaches its 
bulk value. Increasing Dp causes initially a sharp rise in 
AEq. The condensation energy then saturates, exhibit- 
ing damped irregular oscillations, reflecting the competi- 
tion between magnetism and superconductivity. Again, 
in all cases superconductivity favors the AP configura- 
tion. At small I (7 = 0.1), AEq for the P configuration 
vanishes beyond Dp > 4, while in the half-metallic limit, 
AEq is an oscillatory function of Dp. The difference in 
condensation energies between P and AP configurations 
at the same / is a weak function of Dp. 

The irregular oscillatory behavior in Figs. 12131 reflects 
the existence of the characteristic length I — (k^ — fcj) -1 
arising from the difference between Fermi wavectors for 
up and down spins in the F layers. Such oscillatory be- 
havior depends on the relation between dp and I. The 
latter in turn depends on /. Since the ratio If dp depends 
on in a simpler way than on J, the oscillations are 
best studied in Fig. [3J There one can see that at larger 
1=1 (hence smaller £), the characteristic oscillations 
have a shorter spatial period than those at I = 0.1. For 
SFS sandwiches with small Dp one finds 24 oscillations 
in the pair amplitude in F. These can be seen in the left 
edge of Fig. ^ The situation for AEq is much more com- 
plicated, since oscillations in the pair amplitude are only 
indirectly reflected there. 

The strong modifications to the superconducting state 
of the sample should be easily detected in measurements 
of the critical current. These switching effects are also 
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FIG. 4: (Color online) The normalized LDOS (at D F = 10 
and Hb = 0.15) spatially averaged over the S region and 
normalized as explained in the text. Results for both P and 
AP configurations are shown for two I values. 

very well reflected in LDOS results. Thus, we show in 
Fig-E]the LDOS N(z, e) averaged over the S region. The 
results are normalized to the normal state bulk value. 



We display results for two values of /, Dp — 10, and 
Hb = 0.15. One can plainly see the difference between 
P and AP configurations: in the AP case the BCS like 
peaks are much more prominent and the gap fairly well 
defined. In the P case no gap exists, although a weak 
BCS like feature is still visible for / = 1.0, while the 
features flatten out nearly completely at small I (I — 
0.1), when the system is no longer superconducting, as 
seen by the the vanishing of the condensation energy in 
this case, (Fig. 

The enhancement of superconductivity in the AP con- 
figuration is not, as one might naively think, simply a 
consequence of the magnetic polarizations canceling in 
the superconductor. As has been found^i the magnetic 
moment induced in the superconductor by the ferromag- 
netic contacts penetrates into the S material only a few 
Fermi wavelengths. We have verified that this is also the 
case here. Thus, the reasons are more subtle. The weak- 
ening of superconductivity by ferromagnetic contacts de- 
pends in a complicated way on the amplitudes for elec- 
tron scattering (both normal and Andreev) at the inter- 
faces. These are to a greater or lesser strength pair break- 
ing. The penetration depth for Cooper pairs into the F 
material is much smaller than for a normal metal and 
this is reflected in the interface scattering amplitudes. 
Our self consistent calculation shows then, that the su- 
perconducting state (with d$ — £o) can better survive 
proximity to two F contacts that have opposite polariza- 
tions. 
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